Hubbard chains network on corner-sharing tetrahedra: origin of the heavy fermion 

state in LiV204 
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We investigate the Hubbard chains network model defined on corner-sharing tetrahedra (the py- 
rochlore lattice) which is a possible microscopic model for the heavy fermion state of LiV2 04. Based 
upon this model, we can explain transport, magnetic, and thermodynamic properties of LiV204. 
We calculate the spin susceptibility, and the specific heat coefficient, exploiting the Bethe ansatz 
exact solution of the ID Hubbard model and bosonization method. The results are quite consistent 
with experimental observations. We obtain the large specific heat coefficient 7 ~ 222mJ/molK^. 
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I. INTRODUCTION 

Recently, it is discovered that a heavy fermion (HF) 
state is realized in LiV2 04,_a_ transition-metal oxide 
with a cubic spinel structure.El U It has the largest spe- 
cific heat coefficient among d-electron systems, 7 
420mJ/molK^.S The origin of the HF behavior without 
almost localized /-electrons is still an important open is- 
sue. Ab initio band calculations suggest that geometrical 
frustration of the spinel structure in which V sites form 
corner-sharing tetrahedra |Jiejtwork(so-called pyrochlore 
lattice) may be important .LJeI Actually, it is known that 
some geometrically frustrated electron systems in_the. 
vicinity of the Mott transition shows a HF behavior .E3~E3 
However, since the electron density of LiV2 04 is quarter- 
filling, the above idea may not be directly applicable to 
this system. Nevertheless, it is expected that the geo- 
metrical frustration plays some crucial role. On the other 
hand, Anisimov et al. proposed that the t2g band of V 
sites is splitted into localized and conduction parts by 
trigonal crystal field, and thus the system is simulated 
by the Kondo lattice model.tj 

FIG. 1. Some portions of onc-dimensional-like bands 
formed by t2g orbital on the pyrochlore lattice. 

Here, we would like to propose another microscopic 
model for the HF state from a quite different point of 
view. We pay attention to the fact that as was shown by 
band calculations, electrons near the Fermi surface con- 
sist mainly of the t2g orbitals of V ion which is quarter- 
filled, and tha.1_the hybridization with p-electrons of oxy- 
gen is smallE'El We will neglect trigonal field splitting, 
since it is smaller than the band width by ~ 1/10.0 
Then we can see that each t2g orbital on the corners of 
the pyrochlore lattice forms one-dimensional (ID) like 
bands along each edge of tetrahedra (see FIG.l). It 
is expected that the hybridization between these ID 
bands is much suppressed by the geometrical configu- 
ration. In other words, we can say that to reduce the 
geometrical frustration the electronic structure main- 



tains the ID-like character. Actually, the band struc- 
ture obtained from this consideratiorL-Eesembles quali- 
tatively that of ab initio calculations .Era This proposal 
is also supported by the following experimental obser- 
vations, (i) Recent neutron scattering measurements 
shows that there exists antiferromagnetic spiru fluctua- 
tion with a staggered wave vector Q ^ 0.84a*.El This Q 
vector is close to that of the quarter-filled ID Hubbard 
model. (QiD = 2'K/{Adv-v) ~ 0.71a*. dv~v is the 
distance between the nearest neighbor V sites. See ref. 
6) (ii) As temperature is raised, the resistivity increases 
monotonically like ^ T for T > T*, and is not-, satu- 
rated in contrast to /-electron-based HF systems.El Here 
T* is the characteristic temperature which is analogous 
to the Kondo temperature of /-electron-based HF sys- 
tems. The behavior is easily understood, if we identify 
T* with a dimensional crossover temperature. Namely, 
T* is regarded as a crossover temperature from ID to 3D. 
For T > T*, the ID-like electronic structure gives rise 
T-linear resistivityEj (iii) Moreover, the Hall coefficier^ 
measured by Urano et al. is much small for T > T*.u 
Since there exist several Fermi surfaces in this system, it 
is rather difficult to calculate the Hall cocfflcient from a 
microscopic model. However, this experimental fact does 
not contradict with the interpretation that for x > T* 
the system has ID-like character. We show schemati- 
cally the temperature dependence of the Hall coefficient 
and the resistivity of our model in FIG.2. It is noted 
that these properties are not changed qualitatively even 
if the sample is not a single crystal, as in the case of some 
experimental situations. 

Based upon these considerations, we construct the 
network of quarter-filled Hubbard chains in directions 
(1,±1,0), (1,0, ±1), and (0,1, ±1), which is defined on 
the corner-sharing tetrahedra. We assume that the 
chains are coupled weakly with each other at the cor- 
ners of tetrahedra. Because of ID-like structure, elec- 
tron correlation effect is much enhanced, leading to the 
large specific heat coefficient. In this paper, we demon- 
strate this scenario by microscopic calculation. We show 
that the temperature dependence of the spin susceptibil- 
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ity calculated from our model is consistent with experi- 
mental observations, and that the specific heat coefficient 
is enhanced like 7 ^ 222mJ/molK^. A similar idea has 
been recently proposed by Fulde et al. who considered a 
network of ID Heisenberg spin chains with spins 1 and 
f/2.li3 However since LiV204 is metallic, we believe that 
our model is more appropriate to discuss the HF state of 
this system. 
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FIG. 2. A schematic view for the temperature depen- 
dence of the resistivity p and the Hall coefficient Rh of 
the Hubbard chains network model. In high T regions, 
the Luttinger liquid parameter ~ 1. The sign of 
for T < T* depends on the curvature of the 3D Fermi 
surface formed in the low energy scale. 



II. HUBBARD CHAINS NETWORK MODEL 

Based on the consideration presented in the introduc- 
tion, we construct the Hamiltonian of our model in the 
presence of a magnetic field h{xi), 

6 

H = Y. ^iD + E ^"/3cL,c^- - 5MB E SLh{x^), (1) 

a—1 a/3,cr,z i 

where is the Hamiltonian of the ID Hubbard model, 
-f^iD = E cLiCQ<Ti-n + h.c. + ria^.naii- (2) 



i'^ai) annihilation (creation) operator for elec- 
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trons on a chains, and S^^ = cl^^{a;:/2)^^iCacr'f h'>^ 
H^, H^^,...H^ correspond to the chains in the (1, 1,0), 
(0,1,1), (1,0,1), (0,1,-1), (1,0,-1), and (1,-1,0) di- 
rections, respectively. Because of t2g symmetry of d- 
electrons, each chain can not hybridize directly at the 
same points. However, in LiV204, it is possible that the 
hybridization between different chains realizes through p- 
electrons of oxygen. Actually the hybridization between 
the chains should be non-local. This non-locality makes 



the model quite complicated. Thus for simplicity, we 
assume that each Hubbard chains hybridizes at each cor- 
ner of a tetrahedron. The hybridization parameter Vap 
is given by. 



V = Vc 
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(3) 



Then the hybridization gives rise the doubling of ID 
bands generating total 12 bands because of the geometri- 
cal structure of a tetrahedron. In spite of this simplifica- 
tion of the hybridization parameter, the band structure of 
this network model is roughly siniilar to that obtained by 
first principles band calculations I3'H Thus we believe that 
our model captures the important features of IfjV20rt. 
Comparing the band structure with that in refs.0 andcl, 
we choose t = 0.25 eV. We use this parameter in the 
following calculations. 

The issue of dimensional crossovesnfrom ID to 3D has 
been studied extensively so far.ll3~ll3 The basic idea of 
these previous works is to expand the free energy in 
terms of inter-chain couplings, and to apply Landau- 
Ginzburg type arguments. To avoid confusion, we would 
like to stress that our system is not an assembly of par- 
allel chains as was considered in the previous studies, 
but a network composed of chains aligned in six dif- 
ferent directions reflecting the topology of tetrahedra. 
However we can apply the basic method of the dimen- 
sionaLcrossover problem to our model. Following Boies 
et al.ta, we use the Hubbard-Stratonovich transforma- 
tion Vaf^cl^^C/B^^ Vaf3cl„^fpf3cr{Xi) + Vapll}i„{Xi)cp„i + 

ya(}ipaA^i)i'l3a{xi). Averaging over Caai and Caat with 
respect to the action of the ID Hubbard model, we carry 
out the cumulant expansion in Vaf^ and h. Then, we have 
the partition function. 



(4) 



where Zid is the partition function of the ID Hubbard 
model. The effective action is given by. 



S(^l,, ^t) = ^(1) + ^ ^(n) + ^(5-(g)^-(-5))(lD)/j2 
n— 2 q 

5(1) = j dxidX2i^l{xi)[V5{xi-X2) 

-V^G^'-Hxi^x2)]Mx2). (6) 

^("^ = ^E^'""^(^l'^2,...,a;2„)14ftVfea...l^ft„a 
xV';3i<,i(2:i)--V'/32„<T2„(a^2n), (7) 
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wherein = (a;„,t„), and (•••)^^^' is a correlation function 
of the ID Hubbard model. {G(i)(a;)}a^ = G'^\x)Sai3 
with Ga\x) the single-particle Green's function of a 
chain, and c'^a^^xi, ...,X2n) is the connected 2n-particle 



J \(1D) 
■C2n-lC2n)c 



V 



correlation function of a chain, (cic 
is a matrix with elements Vaf3 ■ 

In the presence of the hybridization V, the quasipar- 
ticle weight, which is the hallmark of the Ferim- liq- 
uid state, develops as the temperature is loweredJl3 To 
see this, it is sufficient to consider only up to S^^^ 
term of eq.(||). The pole of the single particle Green's 
function G(/c,£) = G^^'> {k, e)[l - VG^^Hk, e)]'^ deter- 
mines the quasiparticle energy spectrum Ek ~ [2(1 — 
6)vcVs/{vc + Vs)]k = vpk. Here Vc is the velocity 
of holon, and Vg is that of spinon. The quasiparticle 
weight is z ~ 2^/v^{AkVc/ Ep)^ / {vc + Vs), with Afc — 
i4Vv',{vcVs)-'/'E/)y(^-'\ 9 = {K, + l/K,-2)/A, 
the Luttinger liquid parameter of the charge sector, and 
Ep the Fermi energy of the ID Hubbard model. 5*") 
{n > 2) terms are relevant perturbations to this Fermi 
liquid state. As long as spontaneous symmetry breaking 
does not occur, the Fermi liquid state is stable. At suf- 



ficiently high temperatures, we can treat 5'(") (n > 2) 
terms as small perturbations. In the following, we take 
into account only S*^^) term. At low temperatures, effects 
of higher order terms S**^"-* (n > 3) are not negligible. To 
estimate the temperature range in which our approach 
is valid, we apply simple scaling argument to the action 
(B). The scaling equation for the effective coupling of 
S^^\ which we denote (74, is given by, 



dg4 



dln(^; 



= (1 - K,)gi 



(8) 



Then, 174 grows to the order of unity at T ^ 
(y/4t)4/(i-'f^<=)Ti., where Tp - 4<, the band width of 
the ID Hubbard model. Thus our approach is applicable 
for T > To = {V/Uf/^^-^-^Tp. For T < Tq, our ap- 
proximation will be less accurate quantitatively. As we 
will see later, Tq is sufficiently small for the parameters 
we use. 

Using the bosonization rule and the operator product 
expansion for the U{1) Gaussian model and the level-1 
SU{2) Wess-Zumino-Witten modely,we,can show that the 
leading term of S^-^'^ is written as,ll3~Ell 



Here p{q) = E<T,fc cI,fc+gC<,,fc, = Efc 4,fe+<7'^i'*=' 

In a similar manner, we can rewrite the last term of eq.(|^) 

as, 

^ E ^"/5i^/32aEE^"°(9'^)^/3i<Ti,/c-H9<^^ia2V'ft^2,fc/l? 

+ E E^^°('?''^)'^ft-l,fc+9<-2V'/32^2,fcM- (10) 
g~Q OJ 

Here xi^('J' ^ ^^^^ and j}P{q ^ Qo,^) are the uni- 
form and staggered parts of the spin susceptibility, re- 
spectively, for the ID Hubbard model. Qq = -k (2 \vl the 
quarter-filling case. 



III. MAGNETIC PROPERTIES 

In this section we consider the magnetic properties of 
the system. Applying random phase approximation, we 
compute the spin susceptibility, 

= E{x'''(9)(i - r(9)x'^(g))-i}.;3, (11) 

a/3 

where {id^ {(l)} afi = xi°(9a)'^Q/3 with qi = (72 = + 

93 = 94 = % + gz, 95 = 96 = + <?7 = 98 = ^2/ - Iz, 



q9 = qio ^Qx- Qz, and qii ^ qu ^ qx - qy {f (9)}a/3 = 
i^EkG%ik + q)G%{k), with Glp = {VG'^VU. and 

G" ^ [V - t>2G(i)]-i. It is noted that in the de- 
nominator of eq.(O) the single-particle Green's function 
appears. Eq.(|llD implies that as the quasiparticle 
weight develops, the spin-spin correlation between Hub- 
bard chains which is mediated by two particle hopping 
is enhanced. This point is important for the following 
arguments. Although we know the low energy expres- 
sion of G'^'(a;,t), it is_diSicult to carry out the Fourier 
transform analyticallyH3 For simplicity, we approximate 
the anomalous exponent of G''^^{x,t) as 6* w 0. Actu- 
ally 6 is sufficiently small for intermediate strength of U, 
e.g. e = 0.059 {Kc = 0.618) for U/At = 2, and 9 = 0.03 
[Kc = 0.71) for U/At = 1. Then the Fourier transform of 
G(l)(:E,t) is given by. 



G 



-2iA^VcVs ,1 _ . Vsje ~ kv^) 1, 
{2tt)^T{vc-v,) M '2tt{v,-Vs)T' 2' 



^4 2Tr{vc-Vs)T' 2'' 



(12) 



where B{x,y) is the beta function, and ^ is a non- 
universal pre-factor which depends on model parameters. 
Unfortunately we do not know the exact value of A. In 
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the following, to avoid overestimating electron correla- 
tion effects, we use the non-interacting value of A. We 
now calculate the uniform susceptibility xil = 0) from 
eqs.(|ll|) and (|l^). In the numerical calculation, we ne- 
glect the off-diagonal components of & and x, which 
may give subdominat corrections. We choose the pa- 
rameter At = 1 eV, and U = 2 eV, the-.value sug- 
gested from photo-emission measurements. Using the 
Bethe ansatz exact, solution, we can obtain Vc, Vs, and 
;^iD(0) = l/(27ru5).bEa We determine V by fitting the 
uniform spin susceptibility with experimental data. The 
result is shown in FIG.3. We set V/t = 0.4506. For 
this parameter, Tq/Tp ~ 10~^. Thus we can apply 
our results to the wide temperature range including suf- 
ficiently low temperature regions. The spin suscepti- 
bility consists of the Pauli paramagnetic term and the 
temperature-independent Van Vleck term. To fit to ex- 
perimental data, we set the Van Vleck term xvv ~ 
2.1 X 10~'^emu/mol. It is noted that the Curie-like tem- 
perature dependence is not due to the presence of lo- 
cal moment as in the case of dense Kondo systems, but 
caused by the development of 3D-like magnetic correla- 
tion between Hubbard chains. As temperature decreases, 
the single-particle weight becomes large, leading to the 
enhancement of 3D correlation mediated by two particle 
hopping. 



3D Fermi liquid state. 

We calculate the Fermi liquid part 7fl taking into 
account contributions from spin fluctuations. Since in 
the ID Hubbard model the staggered component of the 
spin fluctuation is dominant over the uniform one, the 
antiferromagnetic spin fluctuation gives leading correc- 
tions to the self-energy. This is consistent with re- 
cent neutron scattering measurements. □ The propaga- 
tor for the antiferromagnetic spin fluctuation is ob- 
tained from eq.(|ll|). We calculate numerically the 
Fourier transform of the staggered spin susceptibility 
xV^{x,t) = e*Q«^Im{l/[27r2(x2 - v^^it - iS)'^)^''^^ix^ ^ 
vl{t - iSfY^'^Weit), which is divergent for q = Qq)B 
This divergence is cutoff by the momentum scale at which 
the crossover from ID to 3D system occurs, i.e. q ^ Afc. 
Since the non- universal pre- factor oix^P{x, t) is unknown 
for finite U , we use the pre- factor of non-interacting sys- 
tems to simplify the calculation. xii{q,Lu) has peaks at 
Q = (Qo,0,0), (0,Qo,0) and (Qo/2, Qo/2, Qo/2). Ex- 
panding eq.(lT]) around one of these peaks, we have. 



xii(Q + g,^) 



xoiQ) 



(13) 



Here xo(Q) = {x'°(Q)(l " T{Q)x''' {Q))-']}ii, and 
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FIG. 3. Calculated result of temperature dependence 
of the spin susceptibility (solid line). J?he dots are exper- 
imental data quoted from refs.U andcl. 



IV. SPECIFIC HEAT COEFFICIENT: HEAVY 
FERMION STATE 



Let us consider the specific heat coefficient which char- 
acterizes the HF behavior. According to eq.(^, the free 
energy is given by the sum of the ID part and the Fermi 
liquid part. Thus the specific heat coefficient is expressed 
as, 7 = 71D+7FL- It should be noted that this expression 
is valid only around the crossover temperature T*, since 
in the sufficiently low temperature region the higher order 
terms of S^"'^ which are neglected in our approximation 
will cancel the ID-like contribution, and the conventional 
Fermi liquid result must be recovered. Thus, 7 ~ 7fl in 
the low temperature region where the system is in the 



= = 2[12(feoxi^ + bQXncQxl'' + SbQCoxH'x 

+ 166|coxr (Xs°)' + SOcQxl^'bobQxirxl''] 
xl/(l-Ao), 

= [16(6ox'° + bQxl'')cQxlF + 325QCoxf 

+646|coxi°(x^°)' + 64cQX^°6ofoQxi°X^°] 
xl/(l- Ao), 

^0 = 8(6qx^° + box^f + 166o6Qxl°xf 
+64(6qXs° + boxlF)bobQx}Pxl'', 



(14) 

(15) 
(16) 



with xr = xi"('z - 0), xJ" = x^"(g = Qo), 60 = 

{zV/vs)y{TrvF), bg = {zV/vs)^ln(4t/vFAk)/{TTVF), 
Co = 4:{zV/vsf/{3TT^VF), and CQ = (zV/vs)^ /{l6^TVF'^k^)■ 
l]s^ng eq.(13), we compute the self-energy En and obtain 



the mass enhancement factor 
9Sn 192V2F4xo(Q) 



de 



.^W. Hl + (f +e?h]-4.33. (17) 



Then, we end up with the specific heat coefficient. 



(18) 



This is almost consistent with experimentally observed 
values, i.e. 350 ^ 420mJ/molK^. Since our approx- 
imation underestimates electron correlation effects, the 
above result gives a lower bound for 7. 
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V. SUMMARY AND COMMENTS 

In this paper, we have proposed a microscopic sce- 
nario for the HF state of LiV204. It has been shown 
that magnetic, thermodynamic and transport properties 
of LiV2 04 are well understood in terms of the Hubbard 
chains network model. In the high temperature region, 
the system is in the ID-like state showing anomalous 
transport properties. As temperature is lowered, the 
crossover to 3D Fermi liquid state occurs. The origin 
of heavy fermion mass is ascribed to enhanced electron 
correlation of the ID-like structure. 

The validity of our model completely replies on the 
assumption that the system is on the border of the di- 
mensional crossover. We need more experimental tests 
to confirm the validity of our scenario. Very recently, 
Li-NMR measurement underj-ambient pressure has been 
carried out by Fujiwara et al.EZl They found that the spin 
lattice relaxation rate (TiT)^^ increases, as the applied 
pressure is increased, showing the enhancement of the 
spin fluctuation. This behavior is quite different from 
that expected for usual strongly correlated metals like 
/-electron based HF systems. This experimental fact 
may be easily understood in terms of the dimensional 
crossover phenomena described by our model. 

Finally, we make a brief comment on single-particle 
properties. It is known that in the ID Hubbard model 
at quarter-filling, Vc is much larger than Vs for large U. 
In our case [U /At — 2), Vc = 1.9i and Vg = 0.6t. Thus it 
is expected that for T* <T <^ Ep /k-Q the single-particle 
spectrum have two peaks which correspond to charge and 
spin degrees of freedom, indicating spin-charge separa- 
tion. It may be intriguing to search for spin-charge sep- 
aration behavior in LiV2 04 by photo-emission spectro- 
scope experiment. 
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